Self consistent GW determination of the interaction strength: apphcation to the iron 

arsenide superconductors 
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We introduce a first principles approach to determine tlie strength of the electronic correlations 
based on the fully self consistent GW approximation. The approach provides a seamless interface 
with dynamical mean field theory, and gives good results for well studied correlated materials such 
as NiO. Applied to the recently discovered iron arsenide materials, it accounts for the noticeable 
correlation features observed in optics and photoemission while explaining the absence of visible 
satellites in X-ray absorption experiments and other high energy spectroscopies. 



Many metals, semiconductors, and insulators are well 
described by the "standard model" of solid state physics. 
In this picture the excitations are band electrons, and 
their dispersion can be computed quantitatively in per- 
turbation theory starting from the density functional the- 
ory using the GW method [T] . When this standard model 
fails, we talk about strongly correlated electron systems. 
The presence of strong correlations is debated with each 
new material discovery, as for example in the context of 
the iron pnictide superconductors. On the experimental 
side, controversies arose because optical experiments re- 
vealed significant mass renormalizations WWi while Xray 
absorption, core level spectroscopies and resonant in- 
elastic Xray scattering indicated the absence of satellite 
peaks [SI [SI, which are standard fingerprints of strong cor- 
relations. Photoemission studies indicate that the overall 
bandwidth is narrowed by a factor of two [3 H] but sub- 
stantially larger mass renormalizations are present near 
the Fermi level [9j. Similar controversies arose within 
the first principles approaches to the treatment of cor- 
relations with some theoretical studies supporting the 
notion of weak correlations [TUl - fT^ . while others advo- 
cate a more correlated picture [T5HTS] . To make progress 
on this issues one needs to develop fully ab initio tools 
for addressing the problem of determining the strength 
of correlations and test their predictions against experi- 
ments. 

In this letter we introduce a new first principles 
methodology for evaluating the strength of the corre- 
lations based on the self-consistent GW method. This 
approach has been shown to predict accurate total en- 
ergy [ini HO] , and we expect to obtain reliable estimates 
for the interaction strength since this quantity can be 
thought as a second derivative of the total energy with 
respect to the occupation of the correlated orbitals. We 
test successfully the method on the well studied exam- 
ple of a correlated material NiO, and then we apply it 
to a prototypical iron pnictide BaFe2As2. We find that 
the correlations in iron pnictides are strong, as pointed 
out in Refs. [I5HI8] but unlike earlier studies our ab-initio 
method accounts for the absence of well defined Hubbard 
bands in the spectral functions. Our results are thus in 



excellent agreement with experiment and reconcile the 
results of apparently conflicting spectroscopies. 

We start with the one-particle electron Green's func- 
tion in the solid, G, which is measurable in photoemission 



experiments. We split it into G = G 



0-1 



E, where 



G describes the non-interacting system of electrons, and 
Y, is the frequency dependent self-energy. Both G and 
S are matrices in r, r'. The electrons interact among 
themselves via the Coulomb interaction Vc(r, r') = w~p~\i 
however, the mobile electrons screen it and is therefore 
useful to reformulate the problem in terms of a screened 
Coulomb interaction W defined hy W = Vc/(1 + Veil) 
[H [21] where 11 is the exact polarization function. 

Dynamical Mean Field Theory (DMFT) maps the 
many body problem in the solid to that of a correlated 
atomic shell embedded in an effective medium. The 
medium is described by an energy dependent Weiss field 
Q''\ which obeys the following equation 



g'' 
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Here Giocai a-nd 'Siocai are the local Green's function and 
the local self-energy, respectively. 

The electrons in the renormalized atom feel an effec- 
tive retarded Coulomb interaction U{u!). Just like the 
Weiss field 0° of the atom reflects the delocalizing effect 
of the medium at the single particle level, the Weiss field 
U{uj) captures the screening of the interaction due to the 
presence of the other atoms. The Weiss field at the two 
particle level U{lu) obeys the following relation 



U-' = w-^l, 



H 



local 7 



(2) 



where Tiiocai and Wiocai are the local polarization func- 
tion and the local screened interaction, respectively. The 
bare local propagators G°{oj) and bare interaction C/(w) 
are chosen so as to give the exact Giocai and Wiocai when 
all the local Feynman diagrams are summed up. Eqs. ([l]) 
and (pi) are a version of the extended-DMFT equations 
studied for simplified models in Refs. [221 [23]. The key 
idea of this work is to use this approach to estimate the 
correlation strength in the solid, and illustrate the power 



of the method by a practical reahstic self-consistent im- 
plementation. 

The approach shares ideas with other methods to com- 
pute the local interaction matrix U. Like constrained 
LDA, it defines correlations on a correlated orbital. It 
adopts the philosophy of the constrained RPA method 
[Ml [5S] , which divides the bands into a set that belongs 
to the low energy model, and the rest of the bands, which 
contribute to screening. However, instead of the bands, 
our method uses orbitals to divide the polarization opera- 
tor of the lattice into a local part, involving the correlated 
orbital, and the rest, which screens the local interaction. 

We now describe the steps required for the prac- 
tical implementation of the method and its in- 
terface with LDA-hDMFT [21]: i) We perform 
a fully self consistent GW calculation [T. ii) 

We evaluate Gioc and Wioc using the projector 
P(rr',iLiL2), defined in Ref. EH Giocail.L, = 
/drP(rr',tLiL2)G(rr') and Wiocail,L,:L,L. 
J P{rr,tL4Li)W{rr')P{r'r',tL3L2)drdr'. t is the 
atom index and L — {I, m) is the angular momentum 
index, iii) We evaluate T.ioc{t) = Wiocai{T)Giocai{-T) 
and Uioc{t) = Giocai{T)Giocai{-T). iii) We use Eq. ^ 
to evaluate U{uj), which we now denote by U'-^^ . iv) 
We also evaluate the hybridization function Aqw{i^) 
using Eq. (In) and identity Gqw = i^ — Ei^p — ^gw- 
A.GW contains the coupling of the correlated orbitals 
to the valence states of the system Ai(w), and to 
the semicore states A/f(a;), and can thus be repre- 
sented as A{iuj) — J de[AL{e) + AH{e)]/{iuj — e). In 
LDA+DMFT the hybridization to these semicore state 
is eliminated resulting in A^, which is connected to 
GW hybridization by AQ\y{iuj) = A/^(ia;) — iuja where 
a w / deA/f (e)/£^. This factor is then absorbed by 
rescaling of the field ifi — >■ tfj/^/T+~a and consequently 
the interaction matrix used in the LDA-f DMFT calcu- 
lation becomes Ulda+dmft — U^^ (uj = 0)/(l -I- a)^. 
This renormalization is usually very small, and in 
BaFe2As2 is a « 0.05. 

We use this fully ab initio method to determine the 
interaction matrix strength U'^^ and the occupancy of 
the d orbital n^, which fixes the double-counting correc- 
tion of LDA-f DMFT. With this input, the LDA+DMFT 
method becomes a fully ab-initio method. 

The effective interaction obtained with this method 
is a general symmetric tensor with four indices 

2-,{ini},(7<T' '-^nn,m.3,m2,mxWm4,aV^m3(y'''rm2Cr'Wmia- it IS USC- 

ful to inquire to which extent this interaction can be ap- 
proximated in terms of Slater integrals F^ , where k runs 
over 0, ...21. The optimal determination of this parame- 
ters is done with the projector 



4'' - E 



?7ii,Tn2,m3 ,m4 



1 47r 



HereM.o = (21+1)^ M=2,k=i = 5(2/7)2 andAA,=2,fc=2 = 
(10/21)2. The quality of the projection is excellent 
and can be seen by recomputing the Coulomb repul- 
sion from the Slater integrals and comparing the result- 
ing [/"'o™ with the full U matrix. We mention in pass- 
ing that the naive Hartree-Fock like estimation of Slater 
integrals J = {Umm'mm,')m^m', -^2 = 14/1.625 J and 
F4 — 8.75/1.625 J, can lead to a substantial underes- 
timation of Slater integrals. 

We first test our method in an arc-typical charge trans- 
fer insulator NiO. We get the following static values 
of the Slater integrals Fq = 7.9 eV, F2 = 10 eV, and 
F4 = 6.7 eV. If the Hund's parameter J is computed 
from F2 ( F4 ) we get J(i^2) = 1.16 eV ( J(Fi) = 1.24 eV 
). When these parameter are used in LDA-I-DMFT, the 
agreement between the theory and experiment is very 
good [55]. We mention in passing that when the GW 
screened interaction and polarization are computed from 
the LDA Kohn-Sham states (non self-consistent GW) we 
get slightly smaller interaction strength Fq « 7.2 eV. 



12 
lOH 

8 
6 
4 
2 


0.92 
0,90 
0.88 
0.86 
0.84 
0.82 
0.80, 



-FO 
-F2 

.-Fi 


{^ 


-FQ 
-F2 
-Fi 


y^^^W 


y^ 


■ 


y 




_ — ^ "^ 






(c) ■ 


V (d) 






_^ 


'^.'-'^ 


-J{Fi) 

-J{Fi) 






-{UgW - Uatom)/FO 



20 



10 





0.006 

0.005 

0.004 

0.003 

0.002 

0.001 

0.000 



xK 



GW 



1 (^ima V^k m-i-m-i V^lmi) (>j) 



5 10 15 20 10» IQi 10^ 10^ 



FIG. 1: (Color online) (a),(b) Slater integrals versus Mat- 
subara frequency as computed by fully self-consistent GW 
method, (c) corresponding Hund's coupling strength J, and 
(d) the difference between the GW Coulomb interaction and 
its Slater parameterization. 



Next we turn to the Goulomb repulsion in BaFe2As2. 
Fig. [IJa) and (b) show the frequency dependence of the 
Slater integrals for Fe-Sd orbitals on imaginary frequency 
axis in linear and log scale, respectively. At very high fre- 
quency, the interaction is unscreened and approaches its 
atomic value. The density-density Coulomb interaction 
Fq is strongly screened in the solid, while the higher mul- 
tipoles Fi and F^ are much less energy dependent, and 
almost equal in solid as in the atom. 

The static Coulomb interaction Fq is estimated to be 
no less then 5eV, larger then previously estimated by 
constrained LDA [TT] and constrained RPA [TUj. We 
want to remark that the self-consistency of GW is im- 



properties of the 
LDA+DMFT cal- 
mass renormaliza- 
3 — 5), but also 



portant in this material, because the non-self consis- 
tent version of GW leads to weaker interaction strength 
Fo w3.4eV. 

The higher order multipoles F2 and F4 show only a 
weak frequency dependence. The highest multipole F4 
is less screened then F2, and hence a single number 
J does not parameterize the form of the Hund's cou- 
pling very well, as J(-F'2) 7^ J{Fi) in Fig. [l|c). Fi- 
nally Fig. [ltd) shows that the Slater parametrization of 
the GW Goulomb interaction is remarkably accurate in 
BaFe2As2, with error less then 6%. 

We now turn to the spectral 
BaFe2As2. Earlier 5-band model 
culations [15j displayed important 
tion at low energies (m* jm k, 
showed a sharp lower Hubbard band. We now per- 
form computations by newly implemented charge self- 
consistent LDA-t-DMFT(CTQMC) method, based on 
WIEN2K [Ml, and explained in detail in Ref. [27l The 
GW estimate for the Slater integrals, renormalized by 
1/(1 + a)'^ — 0.91 due to elimination of the hybridization 
with semicore states, are Fq = 4.9 eV, F2 = 6.4 eV and 
F4 = 4.3 eV. We use the standard localized limit double- 
counting, which gives for the valence Ud = 6.2 ± 0.05, 
in perfect agreement with GW estimate Ud = 6.2, hence 
there is no uncertainty in appropriateness of the cho- 
sen double-counting correction. The quasiparticle mass 
renormalization obtained by DMFT \s Z ^ 1/2 for lighter 
x^ — y^ and z^ orbitals and Z ^ 1/3 for heavier a;z, yz 
and zx orbitals. 

In Fig.[2Fa) we plot the spectral function A(k, lu) along 
the path shown in the inset. There are three circular hole 
pockets at F and two electron pockets at M, in agreement 
with ARPES [Sni IS] and LDA. The two smaller pock- 
ets are degenerate and their crossing occurs at 0.1667r/a 
and 0.287r/a, in good agreement with ARPES [9], where 
the pocket size was estimated to 0.14 7r/a and 0.28 7r/a, 
respectively. The Fermi velocities at F point towards X, 
predicted by our method, are ^ 0.45 eVA, more then 
twice smaller then in LDA. The velocity is in reasonable 
agreement with experiment where somewhat different ve- 
locities for the two pockets are estimated to be 0.43 eVA 
and 0.32 eVA 0. We also overlay the ARPES intensity 
from Ref. [S] on our spectral dispersion to emphasize good 
agreement. 

The optical conductivity is also a strong test of the cor- 
relation strength, as pointed out in Ref. [2] The strong re- 
duction of Drude weight and the presence of mid-infrared 
peak at - 0.6 eV was noticed in Refs. O El In Fig. ^h) 
we show optics obtained by LDA and by DMFT, and 
we compare it to experimental results of Refs. [H Al- 
though LDA gives a reasonable order of magnitude for 
optics, it clearly disagrees with experiments in strength 
of the Drude peak (wp « 2.6 eV) and position of the mid- 
infrared peak, coming from the interband transitions. In 
contrary our DMFT results give smaller Drude peak of 



0.6 


(a) 






■71 


0.4 








m 


i^ 


0.2 








p^-w 




0.0 


^ 


^ 


^-/ 


:^ 


P 


-0.2 




^ 






-0.4 




^H 


B 






-0.6 






r^™^ 


^b 






r X M 




r 



7000 
5000 



(b) 



exper. Wp = 1.6eF 
DMFT w, =1.6eV 
LDA td^=2.e,cV 
DMFT-wu„=2.0eV' 




„ 3000 - 
b 

2000 



w(eV) 



FIG. 2; (Color online) (a) A{\i.,uj) of BaFe2As2 at T = 
150 K as computed by LDA-I-DMFT. The inset shows the 
path in momentum space, while the grey inset shows the 
ARPES intensity from Ref. [9!. (b) Optical conductivity of 
LDA-I-DMFT method (DMFT) and its comparison with ex- 
periment of Ref. 3] (exper.). Also shown is the LDA optical 
conductivity and LDA-I-DMFT conductivity with substan- 
tiahy smaller Fq of Ref. [TO] (DMFT-w). The legend con- 
tains the strength of the Drude peak, which is broadened 
only due to electron-electron interactions, (c) Temperature 
dependence of the optical conductivity within LDA-+-DMFT. 



strength ujp = 1.6 eV, in very favorable agreement with 
experiments [31 14] . We also notice similar width of the 
Drude peak in DMFT and experiments, which shows that 
the most important channel for scattering in this mate- 
rial is the electron-electron scattering. Finally, the po- 
sition of the mid-infrared peak, which LDA predicts at 
frequency ^ 1.2 eV, appears around 0.6 eV in DMFT, in 
very favorable agreement with experiment [21 S| ■ 

To demonstrate the sensitivity of the optical conduc- 
tivity to the strength of the correlations, we carried out 
LDA-fDMFT calculation using parameters of Ref. [101 
Fo = 2.69 eV, J = 0.79 eV, and n^ = 6.53. The resulting 



quasiparticle renormalization amplitude is Z ^ 0.6, in 
very good agreement with results of Ref. ,10, In Fig. [2lb) 
we show optical conductivity thus obtained with label 
DMFT-w. We noticed that neither Drude peak weight 
(cjp = 2.0 eV) nor the position of the mid-infrared peak 
(« 1 eV) is in good agreement with experiments, thus 
confirming that BaFe2As2 should not be regarded as 
weakly correlated material. 

Finally, Fig. W[c) shows the temperature dependence 
of optical conductivity as obtained by DMFT. While the 
interband transitions are roughly temperature indepen- 
dent, the Drude peak width and strength is temperature 
dependent, substantially sharpening at 150 K compared 
to 300 K, which is the consequence of the coherence in- 
coherence crossover in this temperature range, discussed 
in Ref. ESI- 
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FIG. 3: (Color online) Total and partial density of states of 
LDA-I-DMFT method compared to LDA density of states. 

Fig. IsFa) show the total density of states (DOS) and 
Fc-3d partial DOS. Panel (b) shows the As-4p partial 
DOS. Comparing LDA-^DMFT partial DOS to LDA par- 
tial DOS, we notice that apart from the renormaliza- 
tion of the low energy quasi-particles, and broadening 
of the high energy features, there is only little differ- 
ence between LDA and LDA-I-DMFT momentum aver- 
aged spectral functions. This is in agreement with X-ray 
absorption spectroscopy [6], where good agreement be- 
tween LDA and the experiments was pointed out. Given 
the strong correlation effects present in optics and low 
energy ARPES, it is unusual that no clear Hubbard-like 
satellites of the atomic like Sd^ state can be identified in 
local density of states. 

The DMFT valence histogram [33], describing the 
probability of finding each Fe-3(i atomic configuration in 
the solid as a function of the renormalized energy of the 
atomic state sheds light on the unusual metallic state of 
the iron pnictides. In a weakly correlated metal, almost 
all the atomic configurations are significantly present in 
the ground state of the solid and their energy vary over 
the scale of the hybridization which represents the band- 
width of the metal. In correlated oxides, on the other 
hand, only a few atomic states in each valence have sub- 



stantial weight, which results in sharp Hubbard bands. 
In BaFe2As2, the probability of the atomic ground state 
with valence A^ = 6, iV = 7 and iV = 5 is only 0.014, 0.01, 
and 0.007, respectively. Other states have smaller prob- 
ability, but remarkably all atomic states with valence 5, 
6, and 7 have finite probability larger then 0.001. The 
large occupancy of the extremely large number of atomic 
configurations is reminiscent of an itinerant system. On 
the other hand, unlike the weakly correlated situation, 
the spread of the multiplets of the N = 5 states, (coming 
from the Slater integrals F2 and F4) is ~ 7eV similarly 
the atomic states with A^ = 6 span an energy range of 
6.5 eV. This scale, represents the width of the Hubbard 
bands and is very large, much larger than the scale of the 
hybridization (2.eV). 

We stress that the absence of clear atomic-like satel- 
lite excitations is not due to weak correlations in FeAs 
materials, as suggested in Refs. fT0Hl2t but rather due to 
the strength of the atomic multiplet splittings and due to 
the broad bandwidth of the highly polarizable As states. 
This situation, arises for the parameters determined from 
the self consistent GW method. It is significantly differ- 
ent from what is found in the oxides, and is captured by 
the charge self-consistent LDA-I-DMFT calculation. 
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